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1. Introduction 


In this work the phrase "wall data". refers to the data acquired from the 
flexible walls of the test section during the normal process of streamlining 
them around a two-dimensional model^. This data comprises the shapes of the 
two flexible walls and their streamwise pressure distributions, together with 
the tunnel reference flow conditions. 

In principle it is possible to derive much steady-state information on the 
model including force components, moment and. shape. The information on shape 
is the aerodynamic shape of the model, that is, its thickness combined with 
the displacement effect of its boundary layer and wake. Experience gained in 
the earliest days of this research (unpublished work in 1974) indicated that 
there was little likelihood of deriving an accurate shape of the model itself 
using wall data when the data carries the usual levels of experimental error. . 

Much higher precision was required. More recent experience in airfoil testing 
suggested that data on wake displacement thickness was available from the walls, 
with the result that effort was put into deriving this data, and any other which 
might be available. 

In the event, it has proved possible to derive lift force and pitching moment, 
in addition to wake thickness. The extraction of lift and pitching moment relies 
on a balance between pressure forces and momentum changes, while after the 
streamlining of the walls their movement apart from their straight positions 
is simply the displacement effect of the wake. 

It will be seen that although at its present stage of development the data 
which the method provides is not all of good quality, its value lies in the fact 
that it provides an alternate source of information on model behaviour, and in 
some testing (for example, the testing of uninstrumented flow visualization 
models) may provide the only source of information. 


2. Principles 

A two-dimensional flowfield is sketched on figure 1 which shows a pair of 
streamlines, one above and one below the airfoil, which will be followed by 
the flexible walls when they are streamlined. Marked in broken lines are the 
upstream and downstream bounds of the test section, arbitrarily defined by 
geometry but assumed far enough from the disturbances created by the model for 
the static pressure to be constant at the reference value. The walls must be 
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stre^lined for this to be true at the exit. These bounds are chosen to 
be vertical. The streamlines with these end boundaries constitute a 
convenient control volume. Pressure forces act all over, but there are 
momentum fluxes across only the entry and exit planes. The force/ 
momentum flux equation for the test gas is 

o P(s)cos6(s)ds + m(Uj^sinai + U 2 sina 2 ) - Fa = 0 (1) 

. 

cv 

where P(s) is the static pressure acting over the boundary, m is the air 
mass flow rate (equal at entry and exit in the absence of any suction or 
blowing in the test section) and Ux, U 2 and ax, 02 are mean velocities 
and flow angles at entry to and exit from the control volume. Fa is the 
force, positive downward, exerted by the airfoil. 

As the test section is long, it is reasonable to assume flow 
perturbations at the entry plane to be small, and the velocity Ux to be 
close in magnitude to U„, the reference velocity for the test section. In 
contrast, no similar assumption is valid in respect of U 2 . A determination 
of U 2 would require a wake traverse, information not available in routine 
testing. Therefore in this analysis it is simply assumed that U 2 equals Ux? 
under most test conditions, and after wall streamlining, the bulk of the flow 
crossing the exit plane has a velocity close to this. The exceptions are the 
wake, and regions affected strongly by shocks. 

Following these assumptions, (1) becomes 


u P(s) cos6(s) ds + mUoo(sinax-sina 2 ) - Fg = 0 (2) 

cv 

The integral is evaluated as the sum of the pressure forces on finite 
intervals of wall length As in the form 

^ Cp.As.cos0 (3) 


where Cp is a local pressure coefficient, and As the wall interval, here 
taken as one inch. Cp and 6 are wall data. The summation is converted 
into a component of lift coefficient by 


Cli 


J Cp.As.cos9 
c 


( 4 ) 


where c is the airfoil chord. 
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The momentum contributions are similarly reduced to the second 
component lift coefficient by 

Ci 2 = ° Ua>(sinai-sina2) (5) 

c 

where = C^i + 0^2 » and is the reference dynamic pressure. 

The pitching moment about the reference point, which is the 
leading edge having a position relative to the test section shown on 
figure 1, is extracted in much the same way. The inlet momenttan vector 
is assumed to lie mid-way between the walls and to have a slope equal to 
the mean wall slope at jack 1. Similarly at the outlet. The wall pressure 
components identified in equation 3, coupled with appropriate moment arms, 
contribute also to pitching moment, as do the additional force components 
Cp.As.sin6. 

When the walls have been streamlined, they are observed to asymtote 
to fixed distances apart, taken to equal the displacement effect of the 
model's wake. The choice of the streamwise station at which this is to 
be measured is rather arbitrary; in the program given later in this report, 
the average wake thickness is derived from jacks 16-18. 

3. The FORTRAN IV Programme 

The listings with identification of variables are given in Appendices 
1 and 2. The code itself contains many comments. 

There are two parts: data input, and data reduction. The data 
input section, generally ISF, contains streamlined-walls data. The 
examples included may be used as program test cases. ISF passes wall 
shapes and pressure coefficients for analysis to the data reduction 
subroutine ISSL. This subroutine can optionally receive wall Mach number 
distributions as input for conversion to pressure coefficients, simply by 
deleting line number 0009. 

ISSL carries out a curve fitting sequence up to line 0165, for 
pressure coefficients and wall displacements. Interpolations are carried 
out to give information on pressure coefficients, wall slopes and displace- 
ments at one-inch intervals. From figure 1 it can be seen that the first 
interpolated data is required at stations 7 and 8 inches. For this 
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purpose the program curve-fits through the data for jacks 1 to 4, giving 
incidentally the values of the wall slopes in line 0004 of subroutine 
RINTER. The same curve-fit is then carried out for jacks 2-5 but the 
interpolation is then required only for the central portion, at 
stations 13 and 14. The latter pattern is continued throughout the test 
section, except at the downstream end where the curve-fit is used not 
only at its centre (stations 39, 40) but also towards the end, at 
stations 42 and 43. Many comments are included in the FORTRAN listing. 

The accumulated data is then used in lines 0166-on to derive lift 
coefficient (real variable RL) , pitching moment (PM) and wake thickness 
(WTH). 

Descriptions of arrays and. variables are given in Appendix 3^ 

4. Program Test Cases 

The first test case is the potential flow data contained in ISFLBW, see 
Appendix 4(a). For this, streamline shapes and pressure distributions have 
been calculated for the streamlines which pass through the wall anchor 
points. The streamlines are formed by combining a free stream and doublet 
with a free vortex centred on the doublet. In order to create the 
displacement effect of a wake, a source is embedded in what would otherwise 
have been the trailing edge of the cylinder formed by the doublet. 

The circulation is chosen such that with a chord of 4 inches the 
lift coefficient is vinity. The run of this program gives Cl = *9964, an 
error of 0.36%. The vortex is located at the quarter-chord point and 
therefore the exact pitching movement coefficient Cm„ is -0.25. ISF 
gives 0*2441, an under-estimate of 2*4%. The wake thickness figure 
represents the displacement effect of the source. 

Two other typical test cases are given, Appendices 4(b) and 4(c). 

These are runs from the transonic self-streamlining wind tunnel, (TSWT) 
on a 4-inch chord NACA 0012-64 section . Both cases are at transonic 
speeds but case (c) has a supercritical top wall and a substantial 
separation on the airfoil upper surfj^e. ' 
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Comparisons Between Airfoil-Derived and Wall-Derived Data in TSWT 

Many streamlined-^alls cases have been analysed, in addition to the 
above two. Comparisons of lift coefficient Cx. are given on Fig. 2 
indicating good agreement. In contrast, the pitching moment coefficient 
data on figure 3 shows a high level of scatter coupled with a general 
over-estimate of magnitude, leading to the conclusion that at tne present 
state of development, this method can only be used as a rough guide. The 
potential-flow test case (Appendix 4(a)) also shows poor agreement between 
the exact and ISF values of pitching moment. 

Wake displacement thickness is also shown on figure 3, with just 
one check point from a wake traverse 9 inches downstream of the airfoil 
trailing edge, roughly inline with jack 16. Only fair agreement is 
indicated, but methods are available and will be incorporated in due 
course for improving the precision of these measurements. 


6. Conclusions 

1. Wall data from a flexible-walled self streamlining test section can 
give good estimates of the lift coefficient of a two-dimensional model. 

2. Estimates of pitching moment and wake thickness need improvement. 
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LISTING OF ISFLBW 


.fort/list:tt: isflbw 

FORTRAN To 002.5 


C 

C 

C 

C 

0001 
i 0002 

! c; 

0003 

0004 

0005 

C 

0006 

0007 

0008 

0009 

0010 

0011 10 
C 
C 
C 
C 

; c 

c 

I c 

! 0012 

c 

c 

c 

0013 


0014 


0015 


0016 


0017 


TWO DIMENSIONAL MODEL t CL 8 WAKE DISPLACEMENT THICKNESS 
FROM WALL DATA. 

WALLS ARE STREAMLINED....... 

TEST CASE LBW. TrB ARE TOF^» BOTTOM WALL CP'S AT JACK POSITIONS. 
DIMENSION Dl<39)»D2<39)»Ul(39>fU2(39)»S<39)»Sl(39)rS2(39) 
DIMENSION X( 19) »T( 19) »B( 19) yU(39) fRMK 19) yRM2( 19) 

CH IS AN ASSUMED CHORD. HT THE TEST SECTION HEIGHT AT ENTRY. 

CH=4. 

CHD=3. 

HT=6. 

A IS ANGLE OF ATTACK k DEGREES r RMO THE REFERENCE MACH NR. 

A=0. 

RM0=0.0001 
Fl=(l + .2*RM0)i!RM0) 

F2=2/(1.4)|<RM0*RM0) 

WRITE(5flO) 

FORMAT(/' TEST CASEJ LIFTING BODY WITH WAKE. INCOMPFtESSIBLE ♦ ' ) 
RM1I-RM2 ARE TOPrBOTTOM WALL DISPLACEMENTS FROM AERODYNAMICALLY 
STRAIGHTrPOSITIVE WHEN WALL ABOOE STRAIGHT. ALL INCHES. ' 

X ARE JACK STREAMWISE STATIONS. TrB TOP r BOTTOM CP'S AT JACKS. 

LINK WITH ISSLrFORLIB 

THE MOST UPSTREAM JACK <NR.l) IS AT STATION X=6 INCHES. 

DATA X/6. r9. rl2. rl5. ria. r20. r21 . r22. r23. r24. r25. r26. r27. r 
C29. r32. r35. r38. r41 . r44./ 

OUR JACK 20 IS NOT INCLUDED BECAUSE IT HAS NO 
STATIC PRESSURE TAPPING. 

DATA T/.01678r .0l907r . 02091 r . 01 778 r -. 01988 r -. 12316 r -. 20562 r 
C~ . 27383 r - . 28947 r - . 25216 r - . 19506 r - . 14482 r - . 1 086 1 r - . 06746 r - . 041 27 r 
C- . 02955 r - . 02303 r - .01887 r - .01599/ 

DATA B/. 03041 r .03939 r .05436 r .08215 r . 13626 r . 16214 r . 1 1957 r . 03852 r 
C.01282r .0434r . 05786 r . 05068 r .03789r . 01 75 r . 0027 r -. 0031 r -. 00545 r 
C-. 00637 r-. 00666/ 

DATA RMl/. 05887 r . 13121 r . 22474 r . 35545 r . 55777 r. 74894 r . 84624 r 
C.92101r .95695r .9531r . 9222 r . 87928 r . 8341 r . 75155 r . 65331 r .S7751r • 
C.51623r .46486r .42064/ 

DATA RM2/.04545r .09748r .15734r .22416r .28063r . 261 01 r . 20584 r 
C. 13301 r .08153r .05533 r .02942 r-. 0052 r~. 04466 r~. 1236 r-. 22558 r~. 30757 r 
C-. 37469 r-. 43108 r-. 47948/ 

CALL SSL < RMO r A r FI y F2 r X r T r B r RM 1 r RM2 y CH r CHD y HT ) 
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APPENDIX 2 

LISTING OF SUBROUTINE ISSL 


,fort/list:tt: issl .. 

FORTRAN IV V02.5 Fri 09-Jul-82 02:35:i9 PAGE 001 


0001 SUBROUTINE SSL<RM0,A»Fl»F2»X»TfB»RMl»RM2»CH»CHD»HT) 

0002 DIMENSION X(19)»T<19)fB(19)»RMK19>»RM2(19) 

0003 DIMENSION D1 <39 ) , D2<39) f U(39> » U1 ( 39 > > U2<39 ) » S (39) » SI (39 ) » S2<39) 
C CHD IS DISTANCE FROM MODEL AXIS OF ROTATION (WHICH IS ON 

C THE TEST SECTION CENTERLINE) TO THE LEADING EDGE. 

0004 AU1=0. . 

0005 A,U2=0. 

0006 RL=0. 

0007 ARM1=0. 

0008 ARM2=0. - , 

0009 GO TO 45 

C BECAUSE FOR THIS TEST CASE THE INPUT WALL DATA IS ALREADY 
C IN CP FORM. 

0010 DO 20 N=l»19 

C CONVERT WALL CENTERLINE MACH'S INTO CP J 

0011 F3=1.+.2*T(N)*T(N) 

0012 T(N)=< < <F1/F3)**3.5)-1. )*F2 

0013 F3=1.+.2*B(N)*B(N) 

0014 B(N)=< < <F1/F3)**3.5)-1.)*F2 

0015 20 CONTINUE 

0016 40 F0RMAT(I3»4F9.5) 

C Y8fX8 ARE DISTANCES OF L.E. ABOVE»AHEAD OF AXIS OF 
C ROTATION OF MODEL. 

0017 45 Y8=CHD)KSIN<A/57.2958) 

0018 X8=CHD*C0S( A/57, 2958) 

C FOUR INTERPOLATIONS FOLLOW : 

0019 DO 250 N4=l,4 

C N4=lr2 FOR TOP, BOTTOM CP INTERPOLATIONS, 

C N4=3,4 FOR TOP, BOTTOM WALL DISPLACEMENT INTERPOLATIONS. 

0020 IF<N4,EQ.l) GO TO 270 

0022 IF<N4.EQ.3) GO TO 200 

0024 IF<N4.EQ,4) GO TO 240 

0026 DO 50 Nl=l,19 

0027 T<N1)=B(N1) 

0028 50 CONTINUE 

0029 GO TO 270 . . .. .. .. 

0030 .200 DO 60 Nl = l,19 

0031 T(N1)=RM1(N1) _ . 

0032 60 CONTINUE 

0033 GO TO 270 ‘ ‘ ' 

0034 240 DO 70 Nl=l,19 

0035 T(N1)=RM2(N1) - .. 

0036 70 CONTINUE 

0037 270 -DO 560 N3=l,4 

0038 X1=X(N3) 

0039 X2=X(N3+1) 

0040 X3=X(N3+2) 

0041 X4=X(N3+3) .. - _ 

0042 Y1=T(N3) 

0043 Y2=T(N3+1) . 

0044 Y3=T(N3+2) 

0045 Y4=T(N3+3) 

0046 CALL C0EFF(X1,X2,X3,X4,Y1,Y2,Y3,Y4,A2VB2VC2)'“" ■ ■ ‘ 

0047 IF(N3.EQ.l) GO TO 410 - 
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FORTRAN 10 

002.5 Fri 09 -Jij 1-;82 02t35:i9 

0049 


IF(N3.EQ.2) 60 TO 450 

0051 


IF(N3.EQ.3) GO TO 490 

0055 


IF(N3.EQ.4) GO TO 530 

0055 

410 

DO 100 IX=6rl2 

0056 


CALL R INTER ( IX ^ X4 r Y4 , A2 » B2 s- C2 s- SA v UA ) 

0057 


S(IX-5)=SA 

0058 


U(IX-5)=UA 

0059 

100 

CONTINUE 

0060 


GO TO 560 

0061 

450 

DO 110 IX==13^15 

0062 


CALL RINTER (IX ^ X4 » Y4 y A2 i- B2 y C2 y SA i* UA ) 

0063 


S(IX-5)=SA 

0064 


U(IX-5)~UA 

0065 

110 

CONTINUE 

0066 


GO TO 560 

0067 

490 

DO 120 IX~16yl8 

0068 


CALL RINTER(IXyX4yY4yA2yB2yC2ySAyUA> 

0069 


S(IX~5)=SA 

0070 


U(IX-5)=UA 

0071 

120 

CONTINUE 

0072 


GO TO 560 

0073 

530 

DO 130 IX=19y2.0 

0074 


CALL RINTER < IX y X4 v Y4 y A2 y B2 y C2 y SA y UA ) 

0075 


S(IX-5)=SA 

0076 


U(IX-5)=UA 

0077 

130 

CONTINUE 

0078 

560 

CONTINUE 

0079 


DO 700 N3==5yl0 

0080 


X1=X(N3) 

0081 


X2=X(N3T1) 

0082 


X3=X(N3+2) 

0083 


X4=X(N3+3) 

0084 


Y1=T(N3) 

0085 


Y2=T(N3+1) 

0086 


Y3=T(N3+2) 

0087 


Y4=T(N3+3). 

0088 


CALL C0EFF(XlyX2yX3yX4y Y1 y Y2y Y3y Y4yA2yB2yC2) 

0089 


X5=N3+16.-X4 

0090 


S ( N3+1 1 ) =3 . >KA2* ( X5>i<X5 ) +2 ♦ *B2*X5+C2 

0091 


U(N3+ll)=T(N3+2) 

0092 

700 

CONTINUE 

0093 


DO 1050 N3=12yl6 

0094 


X1=X(N3) 

0095 


X2-=X(N3+1) 

0096 


X3=X(N3+2) 

0097 


X4=X(N3+3) 

0098 


Y1=T(N3) 

0099 


Y2=T(N3+1) 

0100 


Y3=T(N3+2) 

0101 


Y4=T(N3+3) 

0102 


CALL C0EFF(X1 yX2yX3yX4y Y1 y Y2y Y3y Y4yA2yB2yC2) 

0103 


IF(N3,EQ.12) GO TO 860 

0105 


IF(N3.EQ.13) GO TO 900 

0107 


IF(N3.EQ.14) GO TO 940 


PAGE 002 
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FORTRAN IV 


V02.5 


Fri 09-JU1-82 02J35:i9 


PAGE 003 


0109 


IF<N3.EQ.15) GO TO 980 

0111 


IF(N3.EQ,16) GO TO 1000 

0113 

860 

DO 880 IX=27i-29 

0114 


CALI. RI NTER < I X » X4 » Y4 r A2 » B2 » C2 r SA r UA ) 

0115 


S(IX-5)=SA 

0116 


U(IX-5)==UA 

0117 

880 

CONTINUE 

0118 


GO TO 1050 

0119 

900 

DO 920 IX==30,32 

0120 


CALL R INTER ( I X » X4 y Y4 » A2 B2 » C2 y SA y UA ) 

0121 


S(IX-5)=SA 

0122 


U(IX-5)::=UA 

0123 

920 

CONTINUE 

0124 


GO TO 1050 

0125 

940 

DO 960 IX=33y35 

0126 


CALL RINTER (IX y X4 y Y4 y A2 y B2 y C2 y SA y UA ) 

0127 


S<IX-5)-~SA 

0128 


U(IX-5)--=UA 

0129 

960 

CONTINUE 

0130 


GO TO 1050 

0131 

980 

DO 990 IX=36y38 

0132 


CALL RINTER ( IX yX4yY4yA2yB2yC2ySAyUA) 

0133 


S(IX-5)=SA 

0134 


U(IX-5)=UA 

0135 

990 

CONTINUE 

0136 


GO TO 1050 

0137 

1000 

DO 1020 IX~39y44 

0138 


CALL R I NTER ( IX y X4 y Y4 y A2 y B2 y C2 y SA y UA ) 

0139 


S(IX"5)=SA 

0140 


U<IX-5)=UA 

0141 

1020 

CONTINUE 

0142 

1050 

CONTINUE 


C 

LOAD APPROPRIATE ARRAYS : 

0143 


IF(N4.EQ.2) GO TO 1320 

0145 


IF(N4.EQ.3) 60 TO 1360 

0147 


IF<N4J-:Q,4) go to 1410 

0149 


DO 1300 IP--=ly39 

0150 


U1(IP)=U(IP) 

0151 

1300 

CONTINUE 

0152 

250 

CONTINUE 

0153 

1320 

DO 1340 IPl=ly39 

0154 


U2(IP1)=U(IP1) 

0155 

1340 

CONTINUE 

0156 


GO TO 250 

0157 

1360 

DO 1390 K=ly39 

0158 


S1(K)-S(K)^ 

0159 


D1(K)=U<K) 

0160 

1390 

CONTINUE 

0161 


GO TO 250 

0162 

1410 

DO 1440 K=ly39 

0163 


S2(K)==S(K) 

0164 


D2 ( K ) =U ( K ) 

0165 

1440 

CONTINUE 


C 

INTERPOLATIONS COMPLETE* 
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Fri 09-Jul-82 02:35:19 


FORTRAN IV V02»5 


PAGE 004 


C 

C 

0166 

0167 

0168 

0169 

0170 

0171 

0172 

0173 

0174 

0175 1600 
C 

0176 

0177 

0178 

0179 

0180 
0181 
0182 

C 

0183 

0184 

0185 

0186 

0187 

0188 

0189 

0190 

0191 

0192 

0193 

C 

C 

C 


0194 

1680 

0195 


0196 


0197 


0198 

1705 

0199 


0200 

1710 

0201 


0202 


0204 


0205 

1708 

0206 


0207 

1709 

0208 


0209 


0210 

1698 

0211 

1699 

0212 

1700 


INTEGRATION OF THE CONTRIBUTIONS OF WALL PRESSURE TO LIFT % 
PITCHING moment: 

DO 1600 N=lf39 
A1=ATAN(S1(N) ) 

AU1=AU1+U1(N)*C0S<A1 ) 

A2-=ATAN(S2(N) ) 

AU2=AU2+U2(N)*C0S<A2) 

ARMl=ARMl fUl (N)*C0S(Al)*(24.56“N--5.-XS) 

ARM 1 = ARM 1 --U 1 ( N ) *S IN < A 1 ) ( ( HT/2 ♦ ) fD 1 ( N ) - Y8 ) 

ARM2=ARM2+U2 < N ) *COS ( A2 ) * ( 24 ♦ 56-N-5 » -X8 ) 

ARM2=ARM2fU2<N)*SIN<A2)>ic< (HT/2. )-D2<N)+Y8) 

CONTINUE 

SUM THE COMPONENTS OF LIFT: 

RL=(AU2-AU1)/CH 
SL1=(S1(1)+S2(1) )/2. 

AS1=ATAN(SL1) 

RL=RL+2.*HT*SL1/CH 
SL2=(S1(39)+S2(39) )/2. 

AS2=-ATAN(SL2) 

RL=RL-2.*HT>KSL2/CH 

SUM THE COMPONENTS OF PITCHING MOMENT*. 

Y9=(D2(1)+HT+D1(1) )/2. 

Y1=(18.56-X8)*SL1 
Y2=(Y8f (HT/2. ) )-Y9 
B1 = (Y2-Y1)>!<C0S(AS1) 

Y9=(D2(39)+HT+D1(39) )/2. 

Yl=SL2:j((19.44+X8) 

Y2=(Y8+(HT/2. ) )-Y9 
B2=(Y2fYl)*C0S(AS2) 

BM1=2 . >KHT>KB1/ ( CH*CH ) 

BM2==2 . *HT>KB2/ ( CH*CH ) 

PM=BM2-BM 1 + ( ARM2-ARM 1 ) / ( CH*CH ) 

THE WAKE DISPLACEMENT THICKNESS UTI-I IS TAKEN AS THE AVERAGE 
MOVEMENT APART OF JACKS 16-18: 

WTH= ( RM 1 ( 1 6 ) +RM 1 ( 1 7 ) +RM 1 ( 1 8 ) -RM2 ( 1 6 ) -RM2 ( 1 7 ) -RM2 ( 1 8 ) ) /3 
CALL IDATE(IDlrID2j ID3) 

WRITE (5 » 1700 )RL»PM»UTH 
WRITE(5rl705)ID2rIDl»ID3 

F0RMAT(4Xf "ANALYSIS DATE ' r 313 y 23X r r ' INCHES ' ) 

TYPE 1710 

FORMAT(" WANT INTERPOLATIONS ? INPUT 1 FOR NO:") 

READ(5»>!«)PP 
IF(PP.EQ.l) GO TO 1800 
WRITE(5» 1708) 

FORMAT(" STAT-, DISPLACEMENTS CP SLOPES') 

WRITE(5f 1709) 

FORMAT (' ION UPPER LOWER UPPER LOWER UPPER LOWER 
C' ) 

DO 1699yL=l»39 

WRITE(5yl698)L+5yDl (L) fD2 (L) tUI (L) yU2(L) ySi (L) yS2(L) 
F0RMAT(I3y6F9.4) 

CONTINUE 

F0RMAT(3Xy' CL='yF7.4y' CM='yF7.4y' WAKEDISP. THICKNESS^ ' y F6 . 4 ) 
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FORTRAN IV V02»5 Fri 09-Jul-82 02J35:i9 PAGE 005 

0213 1800 AU2=AU2/CH 

0214 AU1=AU1/CH 

0215 AS1=AS1*2.*HT/CH 

0216 AS2=AS2*2»*HT/CH 

0217 ARMl-=-ARMl/(CH*CH) 

0218 ARM2-ARM2/(CH>!<CH) 

0219 WRITE (5? 1804) 

0220 1804 FORMAT <' BREAKIiOUN t ' 14X ' FLEX* WALL MOMENTUM') 

0221 WRITE(5rl805) 

0222 1805 F0RMAT(18X' PRESSURES FLUXES') 

0223 WRITE(5yl806) 

0224 1806 FORMAT (18X' TOP BOTTOM INLET OUTLET') 

0225 WRITE(5» 1802)“AUlyAU2»ASly-AS2 

0226 WRITE<5y 1810)ARMlyARM2y-BMlyBM2 

0227 1810 FORMAT(' MOMENT COMPONENTS: 'y4F12*4) 

0228 1802 FORMAT</' LIFT COMPONENTS: 'y4F12.4) 

0229 END 


FORTRAN IV V02.5 Fri 09-Jul-82 02:42:i0 PAGE 001 

SUBROUTINE COEFF ( XI y X2 y X3 y X4 y Y1 y Y2 y Y3 y Y4 y A2 y B2 y C2 ) 

COEFFICIENTS A2-C2 FOR CUBIC THRO' THE 4 POINTS ( 0 y 0 ) y ( X5 y Y5 ) y 
♦*(X7yY7) : 

X5-X3-X4 
X6=X2-X4 
X7==X1-X4 
Y5=Y3-Y4 
Y6=Y2-Y4 
Y7=Y1-Y4 

B1:=(X5*X5)-(X5#X5>KX5)/X6 
B3=X7*X7-(X7*X7*X7)/X6 
C1=X5- ( X5X<X5*X5 ) / ( X6*X6 ) 

C3=X7- < X7*X7)KX7 ) / < X6*X6 ) 

Z1=Y5-<Y6*(X5*X5*X5)/(X6*X6*X6) ) 

Z3=Y7-Y6* < X7*X7*X7 ) / < X6*X6>KX6 ) 

C2=( <Z1*B3)/B1)-Z3 
C2=C2/( ( <C1*B3)/B1)-C3) 

B2=(Z1-<C2*C1) )/Bl 

A2=( (Y5-<B2*X5*X5)-C2*X5) )/(X5*X5*X5) 

RETURN 
END 


! FORTRAN IV V02.5 Fri 09-Jul-82 02:43:10 PAGE 001 

0001 SUBROUTINE RINTER ( IX y X4 y Y4 y A2 y B2 y C2 y SA y UA ) 

C VALUES BETWEEN JACKS BY INTERPOLATION WITH CUBIC: 

0002 X5=IX-X4 

0003 UA=Y4+A2*X5*X5*X5+B2*X5>KX5+C2*X5 

0004 SA=(3**A2*X5*X5)+2.>KB2*X5+C2 

0005 RETURN 

0006 END 


0001 

C 

C 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 
0019 


f ^ 



APPENDIX 3 


AREAYS AND VARIABLES USED IN ISF AND ISSL 


PROGRAM :LSFLBU*F-C)R 
ARRAYS} 

X JACK STREAMUISE STATIONS. 

T TOP WALL PRESSURE COEFFICIENTS (OR MACHS). 

B DITTO V BOTTOM. 

RMl TOP WALL MOVEMENTS UP FROM AERODYNAMICALLY STRAIGHT. 

RM2 DITTO V BOTTOM. 

D1 TOP WALL INTERPOLATED DISPLACEMENTS UP. 

D2 DITTO, BOTTOM. 

U;L top wall INTERPOLATED PRESSURE COEFFICIENTS. 

U2 DITTOvBOTTOM. 

S.1 TOP WALL SLOPES AT ONE-INCH INTERVALS. 

S2 DITTO .BOTTOM. 

S,U DUMMY ARRAYS. 

variables: 

A ANGLE OF ATTACK OF AIRFOIL, DEGREES 

CH CHORD 

CHD DISTANCE FORWARD OF LEADING EDGE (WHICH IS 

MOMENT REFERENCE POINT) FROM AXIS 
OF ROTATION(WHICH IS ON THE CENTERLINE). 

F.1. ,F2 COMPRESSIBLE FLOW PARAMETERS. 

HT HEIGHT OF TEST SECTION AT ENTRY. 

RMO REFERENCE MACH NUMBER. 


ALL DIMENSIONS ARE INCHES 



SUBROUTINE 1SSL.,F0R 


ARRAYS : 

T»B THESE CAN INITIALLY CONTAIN WALL MACH NUMBER 

DISTRIBUTIONS»LATER CONCERTED TO PRESSURE 
COEFFICIENT FORM. ARRAY T IS LOADED 
SEQUENTIALLY WITH DATA FOR INTERPOLATION. 

REAL. UAR I ABLEST 

ARMl SUMMATION OF TOP WALL PRESSURE COEFFICIENT 

CONTRIBUTIONS TO PITCHING MOMENT. 

ARM2 DITTO r BOTTOM. LATER IN COEFFICIENT FORM (LINES 

217 f 218). 

AS! FLOW ANGLE AT INLET TO C . U . ( RADI ANS ) . ALPHA 1. 

AS2 DITTO » AT OUTLET. ALPHA 2 ON FIG.l. 

AUl SUM OF CP>lcS*COS(WALL ANGLE). LIFT COMPONENT DUE 

TO PRESSURE ON TOP WALL. S=1 INCH. FORCE 
DOWN IS POSITIVE. 

AU2 DITTO y BOTTOM WALL. FORCE UP IS POSITIVE. 

A1 LOCAL TOP WALL SLOPE AS AN ANGLE (RADIANS). 

A2 DITTO y BOTTOM WALL. 

> THEN A COEFFICIENT IN THE CUBIC EQUATION USED 
IN INTERPOLATION. THE EQUATION HAS THE 
FORM Y-=A2.t ( X ) "3 i B2)i< ( X ) '■'2+C2*X 

BMl CONTRIBUTION TO PITCHING MOMENT COEFFICIENT CM 

OF INLET MOMENTUM FLUX. NOSE UP POSITIVE. 

BM2 DITTO AT OUTLET. NOSE DOWN POSITIVE. 

B1 MOMENT. ARM OF INLET MOMENTUM VECTOR ABOUT MOMENT 

REFERENCE POINT (SEE FIG.l). WHEN 
POSITIVE THE VECTOR PASSES ABOVE THE 
REFERENCE POINT. 

B2 DITTOyBOTTOM. ALSO A COEFFICIENT IN THE CUBIC. 

C2 COEFFICIENT IN CUBIC. 

f;? isentropic compressible flow function. 

PM pitching moment coefficient about leading edge. 

RL LIFT COEFFICIENT CONTRIBUTION FROM WALL PRESSURES. 

SA . SLOPE DERIVED BY CURVE FITTING AND INTERPLOATION . 

SLl ■ MEAN SLOPE OF FLOW AT ENTRY TO C.V. 

SL2 DITTO AT EXIT'. 

UA PRESSURE COEFFICIENT OR WALL DISPLACEMENT DERIVED 

BY CURVE FITTING AND INTERPOLATION. 

X1-X4) DATA TO WHICH THE CUBIC IS TO BE FITTED. 

Y1-Y4) 

X8 DISTANCE OF MOMENT REFERENCE POINT AHEAD OF AXIS 

OF ROTATION. SEE FIG.l. 

Y1 RISE IN THE INLET FLOW VECTOR LINE BETWEEN THE 

INLET % THE MOMENT REFERENCE POINT. LATER 
THE SAME FOR THE OUTLET VECTOR. 

Y8 DISTANCE OF REFERENCE POINT ABOVE THE CENTERLINE. 

Y? HEIGHT OF INLET (OR OUTLET) VECTOR ABOVE THE LOWER 

WALL (WHEN WALL STRAIGHT) WHERE VECTOR 
CROSSES THE INLET (OR OUTLET) BOUNDARY 
OF THE C.V. 


SUBROUTINE COEFF*. 

X5y Y5) 

X6yY6) INPUT DATA REFERENCED TO X4 & Y4. 
X7rY7) 



APPENDIX 4Caj 
TEST CASE 


.RUN ISFLBW 

TEST case: lifting body with wake, incompressible. 

CL= 0.9964 CM--0.2441 WAKE DISP. THICKNESS=0 . £1906 

ANALYSIS DATE 9 7 82 INCHES 

WANT INTERPOLATIONS ? INPUT 1 FOR NO: 


STAT 

- DISPLACEMENTS 


CP 

SLOP 

£S 

ION 

UPPER 

. LOWER 

UPPER 

LOWER 

UPPER 

LOWER 

6 

0.0589 

0.0455 

0.0168 

0.0304 

0.0224 

0.0159 

7 

0.0816 

0.0619 

0.0173' 

0.0332 

0.0232 

0.0169 

8 

0.1055 

0.0792 

0.0181 

0.0361 

0.0247 

0.0178 

9 

0.1312 

0.0975 

0.0191 

0.0394 

0.0268 

0.0187 

10 

0.1592 

0.1166 

0.0200 

0.0434 

0.0294 

0.0196 

11 

0.1902 

0 . 1366 

0.0206 

0.0483 

■ 0.0327 

0.0204 

12 

0 . 2247 

0.1573 

0.0209 

0.0544 

0.0365 

0.0212 

; 13 

0 . 2625 

0.1797 

0.0219 

0.0615 

0.0402 

0.0225 

• 14 ' 

0.3056 

0.2022 

0.0212 

0.0706 

0.0463 

0 . 0223 

15 

0.3555 

0.2242 

0.0178 

0.0821 

0.0536 

0.0215 

16 

0.4133 

0.2493 

0.0150 

0.0997 

0.0622 

0.0233 

17 

0.4803 

0.2694 

0.0039 

0.1184 

0.0720 

0.0163 

18 

0.5578 

0.2806 

-0.0199 

0.1363 

0.0831 

0.0055 

19 

0.6511 

0.2828 

-0.0621 

0.1619 

0.0964 

-0.0082 

20 

0.7489 

0.2610 

-0.1232 

0.1621 

0 . 0984 

-0.0369 

21 

0.8462 

0.2058 

-0.2056 

0.1196 

0.0890 

— 0 . 0656 

'”•*2 

0.9210 

0.1330 

-0.2738 

0.0385 

0.0581 

-0.0687 

23 

0 . 9570 

0.0815 

-0.2895 

0.0128 

0.0162 

-0.0395 

24 

0.9531 

0.0553 

-0.2522 

0.0434 

-0.0195 

-0.0219 

25 

0.9222 

0.0294 

-0.1951 

0.0579 

-0.0394 

-0.0288 

26 

0.8793 

-0.0052 

-0.1448 

0.050'7 

-0.0457 

-0.0377 

27 

0.8341 

-0.0447 

-0.1086 

0.0379 

-0.0440 

-0.0398 

1 28 

0.7914 

-0.0845 

-0.0838 

0.0268 

-0.0413 

-0.0397 

■ 29 

0.7516 

-0.1236 

. -0.0675 

0.0175 

-0.0383 

-0.0384 

30 

0.7157 

-0.1598 

-0.0552 

0.0109 

-0.0342 

-0.0351 

31 

0.6831 

-0,. 1938 

-0.0468 

0.0061 

-0.0312 

-0.0329 

32 

0.6533 

-0.2256 

-0.0413 

0.0027 

-0.0285 

-0.0307 

33 

0.6259 

-0.2549 

-0.0362 

0.0000 

-0.0262 

-0 . 0283 

34 

0.6008 

-0.2821 

-0.0324 

-0.0018 

-0.0242 

-0.0263 

35 

0.5775 

-0.3076 

-0.0295 

-0.0031 

-0.0224 

-0.0246 

36 

0.5557 

-0.3314 

-0.0269 

-0.0042 

“0.0211 

-0.0230 

37 

0 . 5353 

-0.3537 

-0.0248 

-0.0049 

-0.0197 

-0.0216 

38 

0.5162 

-0.3747 

-0.0230 

-0.0055 

-0.0185 

-0.0204 

39 

0.4981 

-0.3945 

-0.0214 

-0.0059 

“0.0176 

-0.0193 

40 

0.4811 

-0.4133 

-0.0201 

-0.0062 

-0.0166 

-0.0183 

41 

0.4649 

-0.4311 

-0.0189 

-0.0064 

-0.0158 

-0.0173 

42 

0.4495 

-0.4480 

-0.0178 

-0.0065 D.0150 

-0.0165 

43 

0.4348 

-0.4641 

"^o.-or69-- 

-<K0066- 

-0.0144 

-0.0157 

1 44 

0.4206 

-0.4795 

-0.0160 

-0.0067 

-0.0139 

-0.0151 


I breakdown: flex, wall momentum 

PRESSURES • FLUXES 

I TOP BOTTOM INLET OUTLET 

j LIFT components: 0.5095 0.3860 0.0574 0.0434 

! MOMENT components: -0.7322 0.4472 0.2626 -0.2216 

14 



APPENDIX 4(b) 


TEST CASE 72 


RUN 1SF72 

RUN 65 ,A-0-A=2.00 DEG. MACH 0.7030 

CL= 0.1799 CM=~0.0093 WAKE DISP. THICKNESS=0 . 0300 

ANALYSIS DATE 000 INCHES 

WANT INTERPOLATIONS ? INPUT 1 FOR NOt 
1 

breakdown: flex, wall momentum 

PRESSURES FLUXES 

TOP BOTTOM INLET OUTLET 

LIFT components: 0.1843 -0.0156 0.0069 0.0042 

MOMENT components: -0.1803 0.1633 -0.0454 0.0531 


. FORT/LIST :TT: 1SF72.F0R 

FORTRAN IV V02.5 Mon 12-Jul-82 00:0i:23 PAGE 001 

i 

C TWO DIMENSIONAL MODEL: CL & WAKE DISPLACEMENT THICKNESS 

C FROM WALL DATA. 

C ....... WALLS ARE STREAMLINED....... 

0001 DIMENSION DK39) ? D2 < 39 ) ^ U1 ( 39 ) » U2 < 39 ) f S < 39 > s- SI < 39 ) fS2(39) 

0002 DIMENSION X ( 19 ) s- T ( 19 ) s- B < 19 ) f U ( 39 ) » RMl ( 1 9 ) y RM2 ( 19 ) 

0003 CH=4. 

0004 CHD==3. 

0005 HT=6. 

0006 RUN=65. 

0007 RM0=.703 

0008 A-2.0 

C RMO IS REF. MACH NR. A IS A-O-A.y DEGREES. . 

0009 Fl==l. + .2*RM0)f<RM0 

0010 F2=2/( 1 .4*RM0*RM0) 

0011 WRITE(5yl0)AvRM0 

0012 10 FORMAT(/' RUN 65 A-0-A= ' y F4 . 2 y ' DEG. MACH ' y F8 . 4 y / > ? 

)!<**^* I 

C RMlyRM2 ARE TOP y BOTTOM WALL DISPLACEMENTS FROM AERODYNAMICALLY 

C STRAIGHT y POSITIVE WHEN WALL ABOVE STRAIGHT. ALL INCHES. 

C X ARE JACK STREAMWISE STATIONS. TyB TOP y BOTTOM MACHS AT JACKS, 

,0013 DATA X/6. y9. yl2. yl5. yl8. y20. y 21 ♦ y 22 . y 23 . y 24 . y 25 , y 26 , y 27 ♦ y 

C29. y32. y35, y38. y41 . y44./ 

0014 DATA T/.7054y .7051y .7071y .7022y .7017y .708y .7215y .7484y .7614y 
C.748y .7248 y .7157 y ,7096 y .7033 y .705 y .7044y .707y . 7096 y .7082/ 

0015 DATA B/.7008y .6988y .7035y.6982y .694y .6908y .6951y .7085y.7178y 

C. 73 y .713 y .7101 y . 6993 ,. 7039 y . 705 y , 7067 y .705 y . 7053 y .7044/ 

0016 DATA RMl/.006y .0147y .0237y .0356y ,0596y .0942y .1269y.1622y.175y 

C. 1639y . 1398 y . 1151 > .0928 y .0648 y . 0385 y . 0268 y . 0245 y . 0181 y ♦ 0108/ 

0017 DATA RM2/.0106y.0149y.0202y.0284y, 0321 y. 0181 y-.0023y-.0285y 

C-. 0445 y -. 0487 y -. 0372y -. 0262y -. 0142 y -, 0071 y-. 0043 y-. 0043 y-. 0091 y 
C-. 0071 y-, 0059/ . 

0018 CALL SSM(RM0yAyFlyF2yXyTyByRMlyRM2yCHyCHDyHT) 
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APPENDIX 4(C) 
TEST CASE 82 


RUN 1SF82 


RUN 172. A-C)-A=2.00 DEG. MACH 0.8480 


CL= 0.0577 CM= 0.0795 WAKE DISP. THICKNESS=0 
ANALYSIS DATE 9 7 82 

WANT INTERPOLATIONS ? INPUT 1 FOR NO t 
1 

breakdown: flex, wall 

PRESSURES 


.0569 

INCHES 


momentum 

FLUXES 


TOP BOTTOM INLET 


OUTLET . 


LIFT components: 
MOMENT components: 


0.0785 -0.0231 

-0.1094 0.1582 


0.0037 -0.0014 

-0.0575 0.0883 


FORT/LIST :TT: 1SF82 

FORTRAN lU 002.5 Fri 09-Jul-82 03:22:57 PAGE 001 


C 

c 

c 

! 0001 
0002 

0003 

0004 

0005 

I 0006 

0007 

0008 

C 

0009 
' 0010 
0011 

0012 10 
I 
C 
C 
C 
C 
C 
C 

0013 

0014 

0015 

0016 
0017 


0018 


TWO DIMENSIONAL MODEL: CL X WAKE DISPLACEMENT THICKNESS 
FROM WALL DATA. 

....... WALLS ARE STREAMLINED ....... 


DIMENSION D1 ( 39 ) y D2 ( 39 ) y U1 < 39 ) U2 < 39 ) y S ( 39 ) y SI < 39 ) y S2 ( 39 ) 
DIMENSION X(19) yT(19) yB(19) yU<39) yRMl(19) yRM2( 19) 

RM0=.848 

A=2. 

CH=4 . 

CHD=3. 


HT=6. 

RUN=172. 

RMO IS REF. MACH NR. A IS A-O-A. 

F1=1.+.2*RM0*RM0 
F2=2/(l .4*RM0*RM0) 

WRITE(5y 10)RUNyAyRM0' 

FORMAT(/' RUN 'yF4.0y' A-0~A= ' y F4 . 2 y ' DEG. MACH ' y F8 . 4 y / ) ? 


RMlyRM2 ARE TOP y BOTTOM WALL DISPLACEMENTS FROM AERODYNAMICALLY 
STRAIGHT y POSITIVE WHEN WALL ABOUE STRAIGHT. ALL INCHES. 


LINK WITH ISSMyFORLIB. 


X ARE JACK STREAMWISE STATIONS. TyB TOP y BOTTOM MACHS AT JACKS. 
DATA X/6. y9. yl2. y 15. y 18. y20. y21 . y 22 . y 23 . y 24 . y 25 . y 26 . y 27 . y 
C29. y32. y35. y38. y41. y44./ 

DATA T/.8498y .8431y .844y .3382y .8311 y .8253y . 8431 y . 899 y . 9901 y 
C.9602y .8779y . 8625 y .8496 y .8401 y . 8475 y . 8374 y . 846 y . 8489y . 8515/ 
DATA B/.8512y .8406y . 8361 y . 8347 y . 8247 y . 8253 y . 8352 y . 8786 y . 9331 y 
C.963y .8957y .8619y .8479y .8413y .8486y .846y . 8452 y . 8449 y . 8489/ 

DATA RMl/.006y .0068y . 01 1 1 y . 0147 y . 0353 y .0733 y .1121 y .1624y .1948y 
C.188y .152 y .1203 y .0929 y .0613 y .0374y . 0316 y . 0305 y . 0238 y . 0169/ 

DATA RM2/.0113y . 0141 y . 01 19 y . 0182 y . 01 18 y -. 0091 y -. 0378 y -. 0767 y 
C-, 1069 y -. 1208 y -. 0992 y -. 0828 y -. 058 y -. 041 1 y-. 0303 y-. 0295 y-. 0285 y 
C-.0268y-.0215/ 

CALL SSM ( RMO y A y F 1 y F2 y X y T y B y RM 1 y RM2 y CH y CHD y HT ) 



Inlet Boundary Moment Reference Outlet Boundary of 



Straight Walls 

Streamlined Walls 


FIG. 1 GEOMETRY OF MODEL AND TEST SECTION 


TWO-DIMENSIONAL TESTING OF NACA 0012-64: COMPARISONS OF 
LIFT COEFFICIENTS DETERMINED FROM WALL AND MODEL DATA. 
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